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INTRODUCTION 


In ultrasonic nondestructive evaluation (NDE) , the nondestructive 
interrogation energy is transmitted by propagating stress waves. The 
parameters which govern the mechanical behavior of materials also 
affect their stress wave propagation characteristics. 

A useful parameter in describing stress wave propagation in a 
medium is attenuation. Attenuation refers to the energy loss associated 
with the decrease in the stress wave amplitude due to both scattering 
and absorption. It has been observed that attenuation can be correlated 
with the tensile strength of flawed glass fiber composites [1] and the 
residual strength of impact-damaged graphite fiber composites [2]. It 
has also been observed that attenuation of as-fabricated graphite fiber 
composites can be correlated with the compression fatigue life [3] and 
flexural fatigue life [4]. In rock structures, attenuation of seismic 
waves is a major topic in exploration seismology [5]. 

The two major sources of ultrasonic attenuation are absorption due 
to viscoelastic effects of the medium and scattering from inhomogeneities 
in the ultrasonic beam path. 

This study considers only attenuation due to scattering by voids 
in the medium. In the theoretical analysis, the geometries of the voids 
considered are spherical and ellipsoidal. The medium is assumed to be 
isotropic and the incident ultrasonic beam is assumed to be a plane 
longitudinal wave of wavelength much larger than the largest dimension 
of the void. 
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THEORETICAL ANALYSIS 


An incident ultrasonic wave is scattered into various directions 
by voids present in the medium. In this analysis, the energy carried in 
the scattered wave is assumed to be entirely lost, thus contributing to 
ultrasonic attenuation. 

It is also assumed that no interaction occurs between neighboring 
voids in the medium containing many randomly distributed voids. Large 
spacings between voids are not required in order to neglect void 
interactions [6,7], It has been observed that longitudinal wave speed 
and attenuation developed using the void non-interaction assumption 
are valid up to 20% of void volume fraction [8], Thus, scattering from an 
individual void can be treated independently for most void volume 
fractions of engineering interest. So, the combined effect of multiple 
voids is simply the arithmetic addition of their individual contributions. 

An obstacle which scatters an incident ultrasonic wave can be 
described by its scattering cross section y [9], The scattering cross 
section of a scatterer is defined as the ratio of the total power 
scattered to the incident power flux. The scattering cross section has 
dimensions corresponding to area. 

The scattering cross sections for a spherical void and an ellipsoidal 
void will be discussed below. Then an analysis is developed for 
evaluating the ultrasonic through-thickness attenuation based on 
scattering. The medium is assumed to be isotropic. The incident wave 
is assumed to be a plane longitudinal wave have a wavelength much larger 
than the largest dimension of the void. 


Scattering from a Spherical Void 


The exact solution for the scattering cross section y s having a 
spherical void of volume V for an incident longitudinal wave with wave 
number k p is [9] 


Y, 



( 1 ) 


where 


8 C 47 T 


j + 40 


2 + 3c 


(4 - 9c ) 


3 2 2 3 9 4 

2-2 - 2 c + 3 c c 


( 2 ) 


and where 



( 3 ) 


- 2 - 



The wave speed of the incident longitudinal wave is Cp, and c s and 
k s are the wave speed and wave number of the corresponding shear wave, 
respectively, having the same frequency as the incident longitudinal 
wave. 

Eqn. (1) gives the exact scattering cross section of a spherical void 
for an incident longitudinal wave. From eqn. (1), it is observed that 
the scattering cross section is proportional to the square of the 
volume of the void and proportional to the fourth power of the incident 
wave number. Because the wave number kp is related to the wavelength 
Xp as 



the scattering cross section is inversely proportional to the fourth 
power of the incident wavelength. 

The scattering cross section y s in eqn. (1) also depends on 
the parameter c which is the ratio of the shear to longitudinal wave 
numbers through the term g c defined in eqn. (2). Fig. 1 shows the 
scattering cross section of a spherical void for incident longitudinal 
waves versus c, based on eqns . (1), (2) and (3). 


Scattering from an Ellipsoidal Void 


No exact solution exists for the scattering cross section of an 
ellipsoidal void. An approximate solution will be derived below for 
the restrictions to be enumerated. 

A schematic illustrating the scattering of ultrasonic wave by an 
ellipsoidal void is shown in Fig. 2. The center of the ellipsoidal 
void coincides with the origin of the x,y,z cartesian coordinate 
system. The ellipsoid can be described by 



where the principal dimensions of the ellipsoidal void are 2a^, 2 a 2 
and 2a^ in the x, y and z directions, respectively. A plane ultrasonic 
longitudinal wave propagating in the z direction impinges upon the 
ellipsoidal void. The incident wave is scattered by the void into 
all directions described by the polar angles 0 and 4> as shown in 
Fig. 2. The medium is assumed to be isotropic. 

The scattering of ultrasonic waves by an ellipsoidal flaw has 
been considered using the so-called Born approximation technique 
[10,11]. The Born approximation represents the first iteration on the 
integral equation formulation of the scattering problem [10]. The 
physical consequence of the use of the Born approximation in the 
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problem considered here is that the true displacement and strain fields 
inside the voids are replaced by the respective fields in the incident 
wave. When the Born approximation is applied to the scattering of an 
incident plane longitudinal wave having wave number kp by a spherical 
void of radius a Q , good agreement is obtained with the exact solution 
given in eqn. (1), in backscattering for wavelengths such that 

k a < 1 . (6) 

p o 

That is, for scattered waves that propagate back towards the incident 
wave and which satisfy eqn, (6) [11] . 

The differential scattering cross section of an ellipsoidal 

void for an incident longitudinal wave based on the Born approximation 
is* [11] 
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where y e is the scattering cross section of an ellipsoidal void for an 
incident longitudinal wave; dft is the differential solid angle which is 
defined by sin 0 d0 d<J> with reference to 0 and <p as shown in Fig. 2 [12] ; 
X and y are the Lame constants of the parent material; and S pp and Sps 
are the shape factors of the differential scattering cross section for 
incident longitudinal wave scattering into longitudinal and shear waves, 
respectively. The shape factors are [11] 


S 

PP 


4tt a^ a£ a^ (sin AK 


PP. 


(AK y 
PP 


AK cos AK ) 
PP EP_ 


S 

ps 


4t r a n a a 0 (sin AK ^ 
1 15 ]>s 

(AK ) 3 

P s 


AK cos AK ) 
ps ps 


( 8 ) 


* Note that the term (4 tt)^ is missing from eqn. (2.12) in [11]. 
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where 


AKpp = kp ^a 2 sin 2 0 cos 2 <}> + sin 2 0 sin 2 <p + a 2 (1 - cos 0) 2 J 
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The scattering cross section y e of the ellipsoidal void can be 
obtained by integrating eqn. (7) over a closed spherical surface centered 
at the flaw as 
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In order to examine the optimal utility of the Born approximation 
as used in [10,11], a further assumption is introduced by restricting 
this analysis to the very long wavelength regime such that 


k a . « 1 (11) 

P i 


where a. is the largest dimension, of the ellipsoidal void. Eqn. (11) 

represents the regime in which the Born approximation is at its best. 

From eqn. (9), it can be observed that AK pp and AK ps are of the order 

of k p a^. Thus, the very long wavelength assumption in eqn. (11) requires 

AK pp and AK ps to be much less than unity. Then, the sine and cosine 

functions in eqn. (8) may be expanded into their respective Taylor 

series where only the first three terms will be retained. 

After the Taylor series expansions are applied to eqns . (8), 

the result is substituted into eqn. (7) , keeping only terms up to 

AK ^ and AK 2 w hich are evaluated using eqn. (9). Then, the resulting 
pp ps’ o-i o 


expression 



is substituted into eqn. (10) and the 


indicated 


integrations are performed to give the scattering cross section of an 
ellipsoidal void y^ as 
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where V is the volume of the ellipsoidal void given by 
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The corresponding shape factors for the scattering cross section 


are 
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where 
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(16) 


Eqn. (12) gives the scattering cross section of an ellipsoidal 
void for an incident longitudinal wave with the very long wavelength 
assumption defined in eqn. (11), and based on the Born approximation. 

The scattering cross section of an ellipsoidal void given in 
eqn. (12) is similar to that of a spherical void given in eqn. (1) 
because both scattering cross sections are proportional to the square 
of the volume of the void and proportional to the fourth power of the 
incident wave number. However, the shape of the void does not 
affect the scattering cross section of a spherical void beyond the 
contribution of its volume, whereas the void shape does affect the 
scattering cross section of an ellipsoidal void through the Spp and 
Sp S shape factors. 

It can be observed from eqns. (12) through (16) that the 
scattering cross section of an ellipsoidal void is the same if the 
dimensions a^ and a 2 are interchanged. This is as expected because 
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interchanging a-^ and a 2 is equivalent to rotating the void about 
the z-axis, which on physical grounds should not affect the scattered 
power of a plane wave propagating in the z direction. 

Additional observations on the effect of the shape of the void 
can be made by considering voids of the same volume V and the same 
dimension a 3 for an incident wave at a fixed wave number k p . Then, 
the only effect of the shape of the void on the scattering cross 
section is through the term f| defined in eqn. (16). For a void with a^ 
equal to a 2 ; that is, a void that presents a circular cross-sectional 
area orthogonal to the incident wave, a specific value of r| is obtained 
from_eqn. (16) and the corresponding y e is obtained from eqn. (12) via 
the S p p and S ps terms defined in eqn, (15) . For a void such that a^ is 
much larger than a2; that is, a void that presents an elongated 
elliptical cross-sectional area orthogonal tq_the incident wave, a. 
larger value of r| and thus smaller values of S pp and S ps and 
correspondingly smaller y e are obtained. Although the two voids have 
the same cross-sectional area (tt ai a2) orthogonal to the incident 
wave, the elongated void has a smaller scattering cross section. This 
is consistent with the expectation that an elongated void will provide 
less attenuation than that of a less elongated void. 

The very long wavel eng th_as sumption in eqn, (11) results in the 
second terms in both S pp and S pg in eqn. (15) being much less than 
the first terms, unity. In fact, when y e is evaluated from eqns. (12) 
through (16) for a range of r| from zero to 1,000,000, for k p a 3 that 
satisfies eqn. (11) and for c equal to 1.5, y e is found to vary only 
by a maximum of 0.4%. 

Thus , for an incident longitudinal wave having a wavelength that 
is very long compared with the largest dimension of the void as stated 
in eqn. (11), the exact shape of the void has little effect on the 
scattered power or the scattering cross section, which depends 
primarily on the volume of the void. In this case, the scattering 
cross section of a spherical void as given in eqn. (1) having the same 
volume as the arbitrarily shaped void can be used. 

One further observation should be made regarding the accuracy 
of the Born approximation. The scattering cross section of an 
ellipsoidal void as given in eqn. (12) and obtained using the Born 
approximation can be used to compute the scattering cross section of 
a spherical void by letting a^, a 2 and a 3 all be equal to a Q , the 
radius of a spherical void. Fig. 3 shows the ratio of the scattering 
cross section of a spherical void using the Born approximation in 
eqn. (12) to that of a spherical void using the exact solution in 
eqn. (1) versus c for two values of kpa Q . The Born approximation 
solution is observed to be insensitive to changes in k p a Q for very- 
long wavelengths satisfying eqn. (11). From Fig. 3, the approximate 
solution obtained using the Born approximation is found to be at most 
56% of that of the exact solution and this occurs at a value of k p a Q 
of zero and a value of c of approximately 1.7. Thus, the accuracy of 
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the Born approximation for evaluating scattering cross sections is 
marginal, even for very long wavelengths. Nevertheless, the observation, 
based on the Born approximation, regarding the negligible effect of the 
exact shape of the void on the scattering cross section for very long 
wavelengths is expected to remain valid. 


Ultrasonic Through-Thickness Attenuation due to Scattering by Voids 


Ultrasonic attenuation can be measured experimentally by the 
through-thickness technique as illustrated schematically in Fig. 4. 

The ultrasonic through-thickness attenuation will be derived in terms 
of the scattering cross section of the voids and the void volume fraction 
of the specimen. 

As shown in Fig. 4, an input ultrasonic transducer is assumed to 
apply a uniform pressure of amplitude p 0 that is sinusoidal in time 
on the surface of a specimen. It is further assumed that the planar 
area of the ultrasonic beam propagates through the specimen of thickness 
Z without alternation, (Such an assumption is obviously an over- 
simplification of the actual situation as a review of [13] will clearly 
indicate. Nevertheless, this model is adequate for the input-output 
relationship to be derived.) Portions of the input ultrasonic energy 
are scattered by the randomly distributed voids, numbering n in the 
ultrasonic beam path. The output transmitted pressure is measured by 
a receiving ultrasonic transducer mounted on the opposite face of the specimen. 

_0i Z 

The output pressure is p Q e where a is the attenuation which is a 
function of frequency. It is assumed that energy carried in the 
scattered waves is lost and that it represents the only energy dissipating 
mechanism considered here. 

The intensity I of a plane longitudinal travelling wave can be 
related to the pressure p as [12] 


I 



(17) 


where p is the mass density of the material and c p is the longitudinal 
wave speed. Thus, the incident ultrasonic intensity I 0 or power flux 
in Fig. 4 is 


2 



(18) 


In accordance with the definition of the scattering cross section 
Y [9], the power lost due to scattering from a single void is yl 0 . 
Assuming no interaction between voids, the total power lost due to 
scattering P Scattered by n voids is 
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( 19 ) 


P ss x E Y 

Scattered o i 


where n is the total number of voids in the ultrasonic beam and where the 
i-th void has scattering cross section 

Substituting eqn. (18) into eng. (19) gives 


Scattered pc 
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Using the planar area A of the ultrasonic beam, the input ultrasonic 
power can be obtained from eqn. (18) as 


In pc 


A . 


( 21 ) 


Similarly, the output ultrasonic power can be obtained as 
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So, the equation for energy balance is 


p — p = p 

In Scattered Out 


(23) 


Substituting eqns. (20) to (22) into eqn. (23) gives 
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which simplifies to 
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Eqn. (25) can be solved for the attenuation as 


a = 




(26) 


Eqn, (26) gives the ultrasonic through-thickness attenuation of 
a specimen based on the scattering cross section y^ of the i-th void, 

1 n \ 

Fig. 5 shows aJl versus the term ( 1 - — Z y. ] based on eqn. (26). As 

1 A i=l V 


indicated earlier, if the excitation frequency corresponds to a very 
long wavelength such that eqn. (11) is satisfied, the scattering cross 
section of each void can be evaluated from eqn. (1) for a spherical 
void having the same volume as the arbitrarily shaped void. 

For the special case when the medium contains n identical voids, 
eqn. (26) can be simplified to 


a = 



(27) 


where y is the scattering cross section of a single void. For a specimen 
containing a void volume fraction V v , the total volume of voids 
encampassed by the volume of the input-output ultrasonic beam is 


nV = V A l 
v 


(28) 


where V is the volume of a single void. Eqn. (28) can be rearranged 
to give the number of voids in the ultrasonic input-output beam as 



AA 
V ■ 


( 29 ) 


Substituting eqn. (29) into eqn. (27) gives 
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(30) 


a = 




Eqn. (30) gives the ultrasonic through-thickness attenuation a 
due to scattering by voids in terms of the void volume fraction for the 
special case when the medium contains randomly distributed identical 

voids. Fig. 5 shows a£ versus the term (i - \ f) based on eqn. (30) . 

It is reemphasized that the attenuation expressed in eqns. (26) and 
(30) is due only to scattering by voids. The total attenuation of the 
medium containing the voids is 


a T = a b + a 


(31) 


where ot T is the total attenuation and a b is the "base” attenuation; 
that is, the attenuation of the material when no voids are present* 
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CONCLUSIONS 


Ultrasonic longitudinal wave through-thickness attenuation of a 
void-containing medium due to scattering has been considered. The 
attenuation was evaluated assuming no interaction between voids. The 
scattered power was assumed to be entirely lost, thus contributing 
to the ultrasonic attenuation. 

The scattered power due to the presence of a void was described 
by the scattering cross section of the void. An exact solution exists 
for the scattering cross section of a spherical void and was given 
in eqn. (1). The scattering cross section of an ellipsoidal void 
was derived in accordance with the Born approximation and was given 
in eqn. (12). A further assumption, restricting the approximate solution 
to very long wavelengths as given in eqn. (11), was made. 

The most striking result, although not totally unexpected, was 
that the exact shape of a void has negligible effect on its scattering 
cross section in the very long wavelength regime as expressed in 
eqn. CH) , The scattering cross sections of spherical and ellipsoidal 
voids are proportional to the square of the volume of the void and 
proportional to the fourth power of the incident wave number. Thus, 
it is suggested that when very long wavelengths are considered, the 
exact solution of the scattering cross section of a spherical void 
of the same volume should be used for arbitrarily shaped voids. 

The ultrasonic through- thickness attenuation was derived in terms 
of the scattering cross section of the individual voids encountered 
by the ultrasonic beam path and was given in eqn. (26) . For the special 
case when a medium contains identical voids, the attenuation was derived 
in terms of the void volume fraction and the scattering cross section 
as given in eqn. (30) . 

This study demonstrates the scattering effects of voids for waves 
of very long wavelengths satisfying eqn. (11) . The results provide 
an analytical basis for the simple modelling of voids in attenuation 
measurements; that is, model arbitrarily shaped voids as spherical voids 
of the same volume. 

Experimental verification of the results are encouraged. Also, 
it is observed that, as implemented here, the Born approximation is of 
somewhat marginal utility in evaluating scattering cross sections of 
voids. Thus, a different approximation procedure [14] may be employed 
to better estimate scattering cross sections, especially for shorter 
wavelengths as expressed by eqn. (6). 
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NORMALIZED SCATTERING CROSS SECTION OF 



1 Scattering cross section of a spherical void for incident longitudinal waves versus ratio of 
shear to longitudinal wave numbers. 



ELLIPSOIDAL 

VOID 


Fig. 2 Schematic of scattering of ultrasonic plane longitudinal wave 
by an ellipsoidal void. 
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RATIO OF SCATTERING CROSS SECTION OF SPHERICAL 
VOID CALCULATED USING THE BORN APPROXIMATION 
TO THAT OF EXACT SOLUTION / X B0RN \ 




ULTRASONIC BEAM INPUT ULTRASONIC 

HAVING PLANAR EXCITATION PRESSURE p 0 



Fig. 4 Schematic illustrating ultrasonic through-thickness testing of 
a specimen containing voids, n of which are in the ultrasonic 
input -output beam path. 
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$ (neper) 
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1 

Fig. 5 Attenuation a of a specimen of thickness £ due to scattering by voids versus ( 1 - — E y ) 

A j „ -i i 

where A is the ultrasonic beam area and Y is the scattering cross section of the i~th void 

y£ 1 

or ( 1 “ V v y ) for the special case when the medium with void volume fraction V contains many 
identical voids in the beam path, each with scattering cross section Y and volume V. 
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